Matterwave localization in disordered cold atom lattices 
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We propose to observe Anderson localization of ultracold atoms in the presence of a random 
potential made of atoms of another species and trapped at the nodes of an optical lattice, with 
a filling factor less than unity. Such systems enable a nearly perfect experimental control of the 
disorder, while the possibility of modelling the scattering potentials by a set of point-like ones allows 
an exact theoretical analysis. This is illustrated by a detailed analysis of the one-dimensional case. 
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Anderson localization Q is an interference phe- 
nomenon occurring in waves propagating in a static dis- 
order: rather than spreading, the wave remains localized 
in a portion of space. For classical waves, locaHzation was 
observed in 2D water waves 0, and in ID Q and 3D Q 
light beams. When the wave corresponds to a quantum 
particle wavefunction, localization in a disordered poten- 
tial (of infinite spatial extension) corresponds to the ex- 
istence of square integrable stationary states at energies 
for which the classical motion is not bounded . 

Indirect evidence of Anderson localization of electrons 
in condensed matter systems (e.g., the conductance de- 
pendence on the temperature) was obtained in 2D and 
3D and also in thin wires, with an interpretation made 
difficult by interaction effects and the presence of a ther- 
mal bath 0. Truly ID condensed matter systems are 
ordinarily subjected to strong interactions, and there is 
presently an active debate about the role of interactions 
in 2D initiated by the observation of a metal-insulator 
transition in 2D electron gases. 

On the other hand, ultracold atomic gases appear as 
very favorable systems for experimental study of Ander- 
son localization of matterwaves. These systems have the 
advantage of being very flexible: due to a very weak 
coupling to the environment they are virtually immune 
to unwanted decoherence while enabling the possibility 
of coupling to a specifically en^eered thermal bath Q 
or an effective magnetic field for the aim of prob- 
ing their effect on the localization. The Feshbach reso- 
nance allows a controllable introduction of interactions 
whose strength can be chosen at will ^'^'^ therefore 
opens the possibility of experimental tests for interaction- 
localization effects as appeared in models 0] which, al- 
though analytically solvable, were considered unrealistic 
until now. Finally, the dimensionality of the gas is also 
adjustable by the use of tailorable optical potentials [l^ . 

A natural way of producing a disordered potential for 
atoms is by using a laser speckle 0, Q . This requires 
however a very careful control of the speckle, to ensure 
that the absence of spreading of a matterwave is not due 
to a trapping of an atom in a local potential minimum. 



which is particularly challenging in 3D where only the 
lowest energy states may be localized. Also, comparison 
with theory generally requires a numerical, rather than 
analytical, solution of Schrodinger's equation. 

In this work, we propose a way to create an almost per- 
fectly well controlled disordered potential for an atomic 
matterwave, which can even be determined by a direct 
measurement. Moreover, this potential can be justifiably 
modelled by point-like scatterers which allows an analyt- 
ical, and often exact [lHll7l| study of the localization and 
its observability. Although we focus on ID, our scheme 
is applicable to any dimension. In fact, in 2D it can be 
subjected to an exact analysis even in the presence 
of a magnetic field. 

We consider a gas of atoms trapped at the nodes of 
an optical lattice of spatial period b (Fig. ^ : each atom 
is cooled down to the ground vibrational state of the lo- 
cal micro-well; each node is occupied by an atom with 
probability p independently of the other nodes, multiple 
occupancies assumed not to occur; tunnelling between 
neighboring sites is made negligible by choosing the mod- 
ulation depth of the optical lattice much larger than the 
energy 1t?{'k /bY /vis, where ttIs is the trapped atom mass, 
to ensure that the spatial configuration is static [l8l |. 

The set of trapped atoms, designated below as 'scat- 
terers', will act as a random potential for atoms of an- 
other species or of the same species but in another inter- 
nal state, denoted as 'test particles'. One should ensure 
that the test particles (unlike the scatterers) will not be 
trapped by the optical lattice. This can be achieved by 
using two different species with sufficiently different res- 
onance frequencies [l9j. In what follows, we assume that 
the test particles experience as an external potential only 
the interaction potential with the scatterers. To ensure 
elastic scattering we also require that the incoming ki- 
netic energy of a test particle is less than the level spacing 
hwg of a trapped scatterer. 
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FIG. 1: Atoms ('scatterers') trapped at the nodes of a peri- 
odic optical lattice and cooled to the vibrational ground state. 
The average occupancy is p (~ 1/3, here). The test particle 
experiences only the delta potential (in black) created by the 
scatterers and is blind to the periodic optical one (in grey). 

test particle, uog is the oscillation frequency of a scatterer 
in the local micro- well of the lattice. Eq.([3l ensures, by 
energy conservation, that at the end of a scattering event 
the scatterer is not left in an excited vibrational level: 
hence, the disordered potential is static. 

Focusing on ID we assume that the test particle is 
strongly trapped in a matterwave guide, with a quantum 
of vibrational energy tuvt much larger than the longitu- 
dinal kinetic energy, so that its transverse y-z motion is 
frozen in the ground vibrational state of the guide. We 
then introduce the model Hamiltonian for the quantum 
motion of a test particle along the lattice direction x: 
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Here the Xj arc the positions of the occupied micro- wells, 
all integer multiple of the lattice spacing b. The effect of 
each scatterer is represented by a Dirac delta potential 
with a coupling constant g. This assumption is reason- 
able when the wavevector k of the test particle is small 
enough: assuming for simplicity that mg = rrit = m, 1/fc 
should be larger than the sizes af\ = {h/mLOs^tY^"^ of the 
harmonic oscillator ground state of a scatterer in a micro- 
well and of the transverse guide ground state of the test 
particle HDl j conditions already ensured by the elasticity 
condition Eq.(^ and by the ID nature of the motion of 
the test particle. A scatterer may then be modelled by a 
zero range potential. When the 3D scattering length, a, 
describing the free space interaction between a scatterer 
and a test particle is much smaller than the harmonic 
oscillator lengths 0^° one is in the so-called Born regime 
and the ID coupling constant g is given by: 
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Our model Hamiltonian was shown in |2l| to lead to 
localization as a consequence of a theorem derived in Ref. 
[2^ . An acceptable quantitative measure of localization, 
that we adopt in this paper, is the decay length of the 
transmission coefficient: taking all the scatterers to be 
in the half space x > and introducing the transmission 
amplitude iAr(A;) of an incoming plane wave of momen- 
tum fc > through a set of scatterers, we define the 



N^+oc Xn — Xi 
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The average (..) is over all possible realizations of the dis- 
order; although, strictly speaking it is not required since 
log |iAr(fc)|, contrarily to ^n(W, itself, is a self-averaging 
quantity for N — > -|-oo |2J, |25j . The transmission and 
reflection amplitudes are given by 



t]si — \/{Rn)ii rN /In ~ ~{R^) 
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where Rn is related to the transfer matrix of the mat- 
terwave through N scatterers, Rn, through the relation 

RN{xi,...,XN) = T{xNy^RNT{xi) (6) 
= GqT{xn ~ XN-i)~^Go . . . T{x2 - xi)'^Go- 

Go is the transfer matrix of a single scatterer at a; = 0: 
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with a = mg/{h?k), and T{x) is the transfer matrix cor- 
responding to a free propagation over an abscissa x: 
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The use oi Rn instead of Rn in Eq. Q does not affect 
\t]^(k)\, but simplifies the calculations, since Rn depends 
only on the variables Xi+i — Xi, which are independent 
random variables with a common probability distribution 
given by: P(Si) = p{l - p)*'"^ where = (xi+i - Xt)/b. 

We calculated the localization constant k numerically, 
by a Monte Carlo averaging over the disorder, taking a 
large enough number of scatterers to ensure convergence 
in Eq.Q). When expressed in units of 1/6, k depends 
on three dimensionless parameters: the filling factor p, 
the reduced momentum kb and the reduced coupling con- 
stant mgb/h^. Assuming for simplicity that ujs = ~ 
so that the harmonic oscillator lengths also coincide. 
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= a^°, and introducing the recoil energy En 



h'^{-Klbf/2m, we find mgb/h^ = 2TT{a/ a^°){huj /2ERf''^ . 
Typically, hw < IOEr. Since we required a/a^° <C 1, one 
should therefore have mgb/fi? ^ 15. In all our numerical 
calculations mgb/fi^ ~ 2.278. 

Figure |21 shows the localization constant dependence 
on the momentum k : Fig. for filling factor p — 0.9 
(solid lines) and p = I (dashed line), and Fig. for 
filling factor p = 0.1 (solid lines). In the case p = 1 
where the scatterers form a finite periodic chain, allowed 
bands (where k = corresponding to infinite localization 
length of the Bloch waves) are separated by forbidden 
gaps. The p = 0.9 case corresponds to a perturbation of 
the periodic chain by an occasional appearance of empty 
sites. Now K takes non-zero values in the former allowed 
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bands while its value in the former forbidden gaps is only 
weakly affected by the disorder. For p = 0.1 the band 
structure of the periodic chain is washed out, except for 
the vanishing of k in the points where kb is an integer 
multiple of tt (corresponding to band edges in the periodic 
case). The peaks and steps on the k curves (marked by 
arrows) will be interpreted analytically below, in terms 
of a phase shift being a rational multiple of tt. 

At a given scattering energy, measuring a finite value of 
K is, strictly speaking, not a proof of localization but may 
be due simply to the presence of a spectral gap as is 
the case in periodic systems. To confirm the existence of 
localized states one should calculate the density of states 
and check that it is also finite. We performed such a 
numerical check, (shown in the upper frames of Figs. [St 
and 133), by imposing periodic boundary conditions in a 
box of size xn — and averaging over disorder. 

We also performed an analytic calculation of k in sev- 
eral limiting cases. The first limiting case is the p —> 
limit for a fixed value of kb and can be treated along the 
lines of Ref. By expanding the matrix product in- 
volving the two factors GoT{xpf — xn^i)~^ and Rn-i in 
Eq.®, one obtains the recursion relation: 



l0g|tAr| = log\totN-l\ 
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where tq and Iq arc the reflection and transmission am- 
plitudes for the transfer matrix Go. When p tends to 
zero, the accumulated phase shift kbs (with P{s) ~ 
p{l — pY~^ as above) between two successive scatterers 
is uniformly distributed between and 2tt (modulo 27r), 
as long as kb is not a rational multiple of tt. Noting that 
;'^| = for \z\ < 1, we obtain 



Cde\og\i 



(l0g|tjv|) = log|to| + (l0g|tAr_i|) 



which leads to 



«;6= (l/2)plog(l + a2). 
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Fig.|2b shows a good agreement of Eq. 1111 with numerics, 
except for the peaks at fc6 = 7r/(s -I- 1) and the steps at 
kb = 27r/(2s-|- 1), for s integer. The calculation is readily 
extended to the p ^ 1 limit, by considering the defects in 
the chain as scatterers on top of a periodic background: 
the propagation of the matter waves in between two con- 
secutive defects of distance sb is given by the transfer 
matrix T{sb) = (T(6)~^Go)'*, and the scattering over a 
defect corresponds to the transfer matrix Gq = T{b)^^. 
Assuming that the incoming energy is in an allowed band 
of the periodic lattice, the matrix T(&)~^Go has eigenval- 
ues of the form e^'^. When 9 is not a rational multiple of 
TT, the procedure of the previous paragraph can be reused, 
replacing Gq by Go and T{sb) by T{sb). We then get a 
formula similar to Ea. (|10|l . with to replaced by ioj which 




FIG. 2: The localization constant k as a function of the in- 
coming wavevector k, for average occupancies (a) p = 0.9, (b) 
p = 0.1. The coupling constant is mgb/f? = 2.278. Black 
solid lines: numerics. Red thin line: analytical result for 
1 — p <C 1 (a), and p <C 1 (b). Peaks and steps in k marked 
by arrows in the inset of (a) and in (b), with the correspond- 
ing values of 6 for (a) and of kb for (b) being rational multiples 
of TT. Blue line in upper frames of (a),(b): density of states. 
Dashed lines in (a): k, for the periodic case p = 1. 

is the transmission amplitude for the transfer matrix Go 
in the basis where T{b) is diagonal. This leads to 



nb ^ {\ — p) log 



exp(ifc6) — p cxp(—ikb) 



l-p 
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where p = | exp(i^) — (1 — ia) exp(ifc6)p/a^ and where 
9 is the solution of cos^ = cos(fc6) + asm{kb) for p < 1. 
This expression agrees well with numerics, see Fig. 
except for the peaks: as shown in this figure, these indeed 
correspond to values of ^/tt that are rationals. 

The third limit we investigated analytically is that of 
a narrow distribution of the phase shift 9 = skb between 
two consecutive scatterers: denoting its average by (0), 
it is assumed that the probability of finding 9 — {9) out 
of an interval of size ^ 1 is small so that kb <1 and that 
the variance = (1 —p){kb/p)'^ < 1. Along the lines 
of [2^ we expand the relation (derived in [2ll|'l: 

Kb^p [ d^pi<P){\og\Mnie)e'^ + Ah2{9)e-'^\) (13) 
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between k and the invariant (Haar) measure ^(0) which 
is a sohition to the Dyson-Schmidt equation: 



(14) 



where M{e) = GoT{ek-^y^. ^{(I),6) is the argument of 
the com plex number Mii(6')e*'^ + Mi2(6')e""^. It is useful 
to note [26j that Eqs. (|13I14|I are invariant under any 9 
independent S'C/(1,1) similarity transformation M{9) — > 
D-^M{9)D. Assuming the energy is in an allowed band 
of a periodic chain of scattercrs of period {9) /k ~ b/p, the 
transfer matrix M{{9)) has then unimodular eigenvalues 
and we choose D so that D-^M{{9))D is diagonal. We 
then expand all functions of 9 in powers of {9 — (9)) up 
to fourth order, which requires an expansion of and 
K to first and second orders in A9^: 

^(^) =^(0)(^) (0) + ... ; k = k(0)+k(1) + «(2) + ... 

(15) 

One has M^°)(0) = l/(27r) and = 0. Denoting ee 
ijg(e»fc6/p(i _ ick)), a, = /m(e'*=^/J'(l - ia)), we get [27| : 



(1)^ _ p(l-p)aV2 ^fc6 
l-a2 



= -p(l-p) 



2 + 



1 - a2 V p ; 



(16) 
(17) 
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In conclusion, we proposed a well controlled way of pro- 
ducing a disordered potential for atomic matter waves, 
by the scattering of test particles on scattercrs trapped 
at the nodes of an optical lattice. We showed how a 
transmission experiment through such a one-dimensional 
disordered chain provides a clear direct evidence of An- 
derson localization at energies where the density of states 
is appreciable. The proposed experiment is one of several 
possibilities, such as the measurement of the (absence of) 
spreading of a wavepacket initially prepared inside the 
disordered medium which is linked to other exactly cal- 
culable 0, 01 aspects of localization: the finite return 
probability and finite inverse participation ratio. The 
proposed scheme is extendable to higher dimensions and 
enables a controlled experimental study of how the lo- 
calization is affected by the introduction of engineered 
thermal-like bath 01 , effective magnetic field , interac- 
tions among test particles jllj . or several coupled chan- 
nels for the transverse motion of the test particle 01 • 
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